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ABSTRACT 


The  development  of  a  reference  test  case  with  repeatable  and  accurate  results  is 
of  paramount  importance  for  the  validation  of  numerical  Fluid  Structure 
Interaction  (FSI)  models  and  simulations.  In  this  study,  a  new  experimental 
facility  was  designed  and  constructed  for  the  investigation  of  FSI  between  fully 
developed  laminar  flow  and  test  section  plate  in  an  enclosed  flow  channel. 
Computational  Fluid  Dynamics  (CFD)  software  was  used  to  model  laminar  flow  in 
a  shallow  rectangular  duct  to  determine  proper  dimensions  of  the  experimental 
apparatus.  The  FSI  problem  was  then  modeled  in  CFD  using  the  predetermined 
dimensions.  The  experimental  set  up  included  construction  of  the  duct,  creation 
of  a  closed  flow  circuit  and  its  instrumentation.  The  metal  test  plate  was  fitted 
with  gages  to  gather  real-time  information  on  the  strain  levels  during  the 
experiment.  Subsequently,  the  experiment  was  performed  and  the  results  were 
compared  with  the  FSI  modeling. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  interaction  between  a  structure  and  a  fluid  in  which  it  is  submerged  or 
a  surrounding  fluid  flow  is  generally  known  as  Fluid  Structure  Interaction  (FSI). 
FSI  problems  exist  in  numerous  applications  of  engineering  such  as  the  stability 
and  response  of  aircraft  wings,  the  flow  of  blood  through  arteries,  the  response  of 
bridges  and  tall  buildings  to  wind,  the  vibration  of  turbine  and  compressor  blades, 
the  oscillation  of  heat  exchangers  [1],  and  underwater  explosion  interaction  with 
submerged  structures  [2-4]. 

With  ever-increasing  computational  power,  there  is  an  ongoing  extensive 
research  on  numerical  methods  for  solving  FSI  problems.  Some  common 
techniques  include  the  finite  element  method  (FEM),  boundary  element  method, 
arbitrary  Lagrangian-Eulerian  (ALE)  and  the  lattice  Boltzmann  method  (LBM). 
These  techniques  were  used  to  solve  fluid  and  structure  equations  either 
individually  or  coupled  [5-8],  and  to  develop  coupling  strategies  for  numerical 
simulations  and  coupling  algorithms  between  computational  fluid  dynamics 
(CFD)  and  computational  structural  dynamics  (CSD)  solvers  [9].  Despite  the 
practical  relevance  of  the  prediction  of  coupled  fluid  and  structural  dynamics, 
much  of  the  research  in  FSI  problems  has  been  numerical  in  nature.  The 
validation  of  numerical  models  for  FSI  simulations  with  reference  to  test  cases 
was  rare. 

In  more  recent  research,  there  has  been  increased  interest  in  the  vibration 
of  plates  either  in  contact  with  or  submerged  in  a  fluid  due  to  the  accessibility  of 
powerful  numerical  methods  in  FSI,  simple  geometry  and  relevancy  in  many 
marine,  geotechnical  and  aerodynamics  applications.  Free  vibration  of  plates 
interacting  with  a  fluid  has  been  investigated  widely  [10-12],  Time-dependent 
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response  of  a  floating  flexible  plate  due  to  an  impulsive  load  [13]  and  high 
frequency  behavior  of  thin  plate  structures  in  contact  with  a  dense  fluid  were  also 
examined  [14]. 

B.  OBJECTIVES 

The  objective  of  this  present  work  is  to  design  and  develop  a  reference 
test  case  for  the  validation  of  numerical  FSI  models  and  simulations,  including 
the  design  and  construction  of  experimental  apparatus.  The  interaction  between 
a  simple  plate  and  a  fluid  flow  under  laminar,  well-defined,  uniform  and  fully 
developed  inlet  and  outlet  flow  conditions  was  investigated.  Phases  of  this  work 
included  modeling  of  fully  developed  laminar  flow  through  a  duct  using  CFD 
software,  determining  the  appropriate  dimensions  of  the  duct  and  velocity  of  flow, 
construction  of  the  experimental  apparatus,  forming  a  closed  flow  circuit, 
measurement  and  collection  of  strain  levels.  The  findings  of  this  work  will  serve 
as  a  baseline  test  case  and  database,  which  can  be  easily  used  to  verify  and 
validate  different  numerical  and  computational  FSI  models.  A  simple  plate  and 
duct  geometry  was  used  intentionally  to  enhance  the  reproducibility  of  the  results 
and  the  measurement  of  moderate  structure  motion. 
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II.  DESIGN  OF  EXPERIMENTAL  APPARATUS 


A.  COMPUTATIONAL  MODEL 

As  an  initial  step  in  the  design  of  the  experimental  apparatus,  fluid  flow 
through  a  rectangular  duct  was  modeled  using  the  finite  element  analysis 
simulation  code,  ANSYS  13.0.  The  main  objective  of  these  simulations  was  to 
determine  the  location  at  which  fully  developed  flow  would  be  achieved,  hence 
determining  the  appropriate  dimensions  for  the  actual  experimental  apparatus. 
This  was  done  through  the  examination  of  flow  velocity  profiles  at  different 
planes. 

The  dimensions  used  for  the  duct  were  a  length  of  2  m,  width  of  0.2  m 
and  height  of  0.05  m,  see  Figure  1. 


Figure  1 .  Rectangular  Duct  for  Initial  Simulation 


As  the  requirement  of  the  project  was  in  the  laminar  regime,  a  fluid  of  high 
dynamic  viscosity  was  chosen  to  attain  Reynolds  numbers  (Re)  of  less  than 
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2300.  Propylene  Glycol  with  a  density  of  965.3  kg/m3  and  dynamic  viscosity  of 
0.06  kg/m  s  was  used  in  the  simulations.  A  uniform  inlet  velocity  ranging  from  0.3 
to  0.8  m/s  was  specified  as  the  inlet  boundary  condition  in  the  +X  direction  for 
different  runs  of  the  simulation,  while  an  outlet  boundary  condition  of  0  Pa  was 
set  as  the  reference  pressure.  The  flow  velocity  profiles  at  multiple  points  of  the 
duct  were  taken  from  three  planes,  XZ  plane,  XY  plane  at  center  and  XY  plane  at 
side  as  shown  in  Figure  2. 


(continued) 
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Figure  2.  Planes  at  Which  Velocity  Profiles  Were  Taken 


For  illustration  purposes,  results  taken  for  one  of  the  runs  (uniform  inlet 
velocity  of  0.8  m/s)  are  shown  in  Figures  3  to  8. 
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Flow  Velocity  Profile  at  X=0.6m 
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Figure  3.  Flow  Velocity  Profiles  at  XZ  Plane  of  Rectangular  Duct 
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Flow  Velocity  Profile  at  X=!0m 


Figure  4.  Flow  Velocity  Profiles  at  XY  Plane  at  Center 

of  Rectangular  Duct 
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Figure  5.  Flow  Velocity  Profiles  at  XY  Plane  at  Side 

of  Rectangular  Duct 
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Figure  6.  Velocity  Vectors  at  XZ  Plane  of  Rectangular  Duct 
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Figure  7.  Velocity  Vectors  at  XY  Plane  Center  of  Rectangular  Duct 
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Figure  8.  Velocity  Vectors  at  XY  Plane  at  Side  of  Rectangular  Duct 

As  shown  in  the  preceding  group  of  figures,  the  flow  in  the  duct  was  able 
to  achieve  near  fully  developed  flow  in  most  parts  of  the  duct  except  at  Z=0.2m, 
with  a  velocity  magnitude  variance  of  less  than  2%. 

A  refined  duct  with  length  of  2.95  m,  width  of  0.25  m  and  height  of  0.032  m 
was  modeled  in  the  next  set  of  simulations.  Diverging  and  converging  sections 
with  length  of  0.15  m  and  end  nose  diameter  of  0.025  m  were  added  at  the 
entrance  and  exit  respectively  to  represent  pipe  connections  from  a  pump  and  to 
a  reservoir,  see  Figure  9. 
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Figure  9.  Duct  with  Diverging  and  Converging  Sections 

The  same  fluid  and  range  of  uniform  inlet  velocities  were  used  for  the 
simulations.  The  flow  velocity  profiles  were  taken  at  multiple  points  at  the  planes 
shown  in  Figure  2,  with  more  points  taken  before  and  after  the  half  of  the  duct. 
The  results  obtained  are  shown  in  Figures  10  to  15. 
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Flow  Velocity  Profile  at  X=1.37m 


Flow  Velocity  Profile  at  X=1.47m 


Flow  Velocity  Profile  at  X=1 .42m 


Flow  Velocity  Profile  atX=1i2m 


(continued) 
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Flow  Velocity  Profile  at  X=1.57m 
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Figure  10.  Flow  Velocity  Profiles  at  XZ  Plane  of  Duct  with 
Diverging  and  Converging  Sections 
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Flow  Velocity  Profile  at  X=0.37m 


Flow  Velocity  Profile  at  X=1.47m 
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(continued) 
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Flow  Velocity  Profite  at  X=157m  Row  Velocity  Profile  at  X=1.62m 


Flow  Velocity  Profite  at  X-1 14m  Flow  Velocity  Profite  at  Xs2.21m 


Flow  Velocity  Profite  at  X-2£8m 


Figure  1 1 .  Flow  Velocity  Profiles  at  XY  Plane  at  Center  of  Duct  with 

Diverging  and  Converging  Sections 


16 


Flow  Velocity  Profile  at  X=0.37m 


Flow  Velocity  Profile  at  X=0.74m 


Flow  Velocity  Profile  at  Xs1.10m 


Flow  Velocity  Profile  at  Xs1.32m 


Flow  Velocity  Profile  at  Xs1.37m 


Flow  Velocity  Profile  at  X=1.47m 


Flow  Velocity  Profile  at  X=1.42m 


Flow  Velocity  Profile  at  X=1.52m 


(continued) 
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Figure  12.  Flow  Velocity  Profiles  at  XY  Plane  at  Side  of  Duct  with 

Diverging  and  Converging  Sections 
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Figure  1 3.  Velocity  Vectors  at  XZ  Plane  of  Duct  with  Diverging 

and  Converging  Sections 


Figure  14.  Velocity  Vectors  at  XY  Plane  Center  of  Duct  with  Diverging 

and  Converging  Sections 
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Figure  1 5.  Velocity  Vectors  at  XY  Plane  at  Side  of  Duct  with 
Diverging  and  Converging  Sections 


It  can  be  seen  from  the  various  flow  velocity  profiles  that  the  flow  comes  to 
near  fully  developed  flow  as  it  approaches  the  halfway  point  along  of  the  duct.  On 
the  other  hand,  from  the  velocity  vectors  plots  as  shown  in  Figure  13,  it  is  clear 
that  there  is  a  concentration  of  flow  through  the  center  of  the  duct  in  the  initial 
portion  of  the  duct  as  well  as  towards  the  end.  This  is  due  to  the  added  diverging 
and  converging  inlet  and  outlet. 

In  the  next  iteration  of  design  and  simulation,  the  objective  was  to  achieve 
a  further  distributed  flow  as  early  as  possible.  In  an  attempt  to  distribute  the  flow 
from  the  center,  two  sets  of  metallic  honeycomb  structures  were  added  into  the 
flow  channel.  The  entrance  portion  was  split  into  a  3-pipe  inlet  instead  of  one 
diverging  section,  as  seen  in  Figure  16. 
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Figure  16.  Honeycomb  Structures  and  Split  Pipe  Entrance 

The  same  fluid  was  used  for  all  of  the  simulations.  The  uniform  inlet 
velocity  at  each  pipe  inlet  was  set  as  1  m/s.  The  flow  velocity  profiles  were  taken 
at  multiple  points  along  the  planes  shown  in  Figure  2,  with  more  points  examined 
before  and  after  the  midpoint  of  the  duct.  The  results  obtained  are  shown  in 
Figures  17  to  22. 
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Flow  Velocity  Profile  at  X=0.74m 
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(continued) 
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Flow  Velocity  Profile  at  X=1 .57m  Flow  Velocity  Profile  at  X*1 .62m 
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Flow  Velocity  Profile  at  X=1.84m 
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Figure  17.  Flow  Velocity  Profiles  at  XZ  Plane  of  Duct  with 

Honeycomb  Structures 
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Flow  Velocity  Profile  at  X=0.37m 


Flow  Velocity  Profile  at  XK)J4m 


Flow  Velocity  Profile  at  X=1.1  Om  Flow  Velocity  Profile  at  X*1.32m 


Flow  Velocity  Profile  at  X=U7m  Flow  Velocity  Profile  at  X=1 .42m 


(continued) 
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Flow  Velocity  Profile  at  X=2<58m 


Figure  18.  Flow  Velocity  Profiles  at  XY  Plane  at  Center  of  Duct  with 

Honeycomb  Structures 
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Flow  Velocity  Profile  atX-!84m 
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Figure  1 9.  Flow  Velocity  Profiles  at  XY  Plane  at  Side  of  Duct  with 

Honeycomb  Structures 
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Figure  20.  Velocity  Vectors  at  XZ  Plane  of  Duct  with  Honeycomb  Structures 
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Figure  21 .  Velocity  Vectors  at  XY  Plane  at  Center  of  Duct  with  Honeycomb 

Structures 
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Figure  22.  Velocity  Vectors  at  XY  Plane  at  Side  of  Duct  with 

Honeycomb  Structures 


From  the  velocity  vectors  plot  in  Figure  20  and  other  velocity  profiles,  it 
was  observed  that  the  flow  did  not  concentrate  at  the  center  and  was  near  fully 
developed  flow  as  it  approaches  the  halfway  point  of  the  duct.  This  showed  that 
the  change  of  diverging  entrance  into  a  3-pipe  inlet  and  addition  of  honeycomb 
structures  were  effective  in  distributing  the  flow  further. 

In  the  next  iteration  of  design  and  simulation,  the  overall  dimensions  and 
honeycomb  structures  were  retained.  Changes  were  made  to  the  entrance  to 
further  split  it  into  a  4-pipe  inlet  and  the  end  converging  section  was  modified  into 
a  drainage  point  in  order  to  simplify  the  final  design.  Both  inlet  and  outlet  sections 
are  0.1524  m  in  length,  giving  a  total  duct  length  of  3.05  m,  see  Figure  23. 
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Figure  23.  Final  Duct  Design 

The  same  set  of  parameters  and  conditions  were  used  to  run  the 
simulations  and  the  results  obtained  are  shown  in  Figure  24  to  29. 


30 


Flow  Velocity,  u(mfs) 


Flow  Velocity,  u  (mis) 


Flow  Velocity,  ufmi'sf 


Flow  Velocity  Profile  at  X=0J4m 


IMS 

0.1 

DOS 

I  0 

N 

-0  05 
•01 
-015 


1  02  04  06 

08/ 

1 

■ - 

015 

01 

0.05 

I  0 

M 

•0  05 
•0.1 
■CIS 


015 

01 

0.05 

I  0 

N 

-005 

•01 

-015 


Flow  Velocity,  u  (m Is) 


015 

01 

005 

1  0 
N 

•0  05 
•01 
-015 


Flow  Velocity,  u  (m/s) 


Flow  Velocity  Profile  at  X=!32m 


0.2 


*04 


oe 


Flow  Velocity,  u  |m/s) 


Flow  Velocity  Profile  at  X=1.42m 


Flow  Velocity.  u(n»/s) 


Flow  Velocity  Profile  at  X=1.52m 


_ 

~~  

L 

02  04  0  6  0. 

1  1 

1 

— — • 

_ 

\ 

■ 

)  02  04  06  0 

1 

X 

— 

Flow  Velocity,  u  (m/s) 


(continued) 
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Flow  Velocity  Profile  at  X=1.57m 
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Figure  24.  Flow  Velocity  Profiles  at  XZ  Plane  of  Final  Duct  Design 
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Flow  Velocity  Profile  at  X=0.37m 
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Flow  Velocity  Profile  at  X^liTm 
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Figure  25.  Flow  Velocity  Profiles  at  XY  Plane  at  Center  of  Final  Duct  Design 
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Flow  Velocity  Profile  at  X=0.37m 
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Flow  Velocity  Profile  at  X=157m 
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Figure  26.  Flow  Velocity  Profiles  at  XY  Plane  at  Side  of  Final  Duct  Design 
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Figure  27.  Velocity  Vectors  at  XZ  Plane  of  Final  Duct  Design 


Figure  28.  Velocity  Vectors  at  XY  Plane  at  Center  of  Final  Duct  Design 
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Figure  29.  Velocity  Vectors  at  XY  Plane  at  Side  of  Final  Duct  Design 

From  the  velocity  vectors  plot  in  Figure  27  and  various  velocity  profiles,  it 
was  observed  that  the  flow  was  well  distributed  and  was  near  fully  developed 
flow  as  it  approaches  the  half  of  the  duct. 

B.  SELECTION  OF  MATERIALS 
1.  Duct 

To  allow  visualization  of  the  flow,  transparent  Polymethyl  methacrylate 
(PMMA,  Plexiglass)  of  thickness  0.0127  m  was  chosen  to  build  the  duct,  as  seen 
in  Figure  30. 
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Figure  30.  Plexiglass  for  Construction  of  Duct 


Plexiglass  has  other  advantages  such  as  ease  of  cutting,  joining,  low  cost 
and  it  is  inert  to  the  chosen  test  fluid  as  well.  Due  to  its  transparency,  the  duct 
could  be  used  for  future  work  such  as  flow  investigation  using  laser  optical 
measurement  techniques. 

2.  Test  Plate 

As  the  experiment  would  be  a  low-pressure  application,  Aluminum  alloy 
(Al  1100)  with  a  minimum  99%  aluminum,  was  chosen  to  be  the  test  plate 
material.  It  has  a  Young’s  Modulus  of  68.9  GPa,  Poisson  ratio  of  0.33  and 
density  of  2710  kg/m3,  as  shown  in  Figure  31.  The  test  plate  was  0.254  m  by 
0.305  m  and  0.000508  m  thick. 


Figure31.  Al  1 100  Test  Plate 
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3. 


Test  Fluid 


As  previously  mentioned,  the  requirement  of  the  project  was  in  the  laminar 
regime,  a  fluid  of  high  dynamic  viscosity  needs  to  be  chosen  as  the  test  fluid  to 
attain  Reynolds  numbers  (Re)  of  less  than  2300.  In  this  case,  Propylene  Glycol 
with  a  density  of  965.3  kg/m3  and  dynamic  viscosity  of  0.06  kg/m  s  was  selected 
as  the  test  fluid. 

4.  Piping  System 

For  flexibility  and  simplicity,  Polyvinyl  chloride  (PVC)  pipes  and  connectors 
were  chosen  to  build  the  piping  system,  (see  Figure  32).  In  addition,  PVC  is 
lightweight,  durable,  low  cost  and  has  low  reactivity.  The  pipes  and  connectors 
were  assembled  using  solvent  welding. 


Figure  32.  PVC  Pipes  and  Connectors 


C.  APPARATUS  SETUP 

1.  Construction  of  Duct 

The  entire  duct  consists  of  three  center  sections  that  were  each  0.914  m 
long,  0.254  m  wide  and  0.0318  m  high,  and  two  sections  for  the  entrance  and 
exit  sections  that  were  each  0.152  m  long,  0.254  m  wide  and  0.0318  m  high. 
Each  section  was  made  out  of  0.0127  m  thick  Plexiglass  plates  assembled  by 
solvent  welding.  For  examples,  see  Figure  33. 
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(a)  Center  Section  of  Duct 


(b)  Entrance  Section  of  Duct 


Figure  33.  Sections  of  Duct 

As  previously  mentioned,  to  further  distribute  the  flow,  two  sets  of 
honeycomb  structures  were  placed  0.152  m  from  the  entrance  and  exit 
respectively.  The  honeycomb  structures  were  made  of  aluminum  and  were  0.254 
m  wide,  0.0508  long  and  0.0318  high,  see  Figure  34. 


Figure  34.  Honeycomb  Structures 

To  ensure  that  the  test  plate  lied  in  the  region  where  the  flow  before  and 
after  the  test  plate  was  fully  developed,  it  was  placed  at  1.57  m  downstream  from 
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the  entrance  of  the  duct.  A  0.203  m  by  0.203  m  window  section  was  cut  from  the 
bottom  plate  of  the  duct  and  the  test  plate  was  joined  to  the  Plexiglass  plate, 
covering  the  section,  using  3M™  super  77  Multipurpose  Spray  Adhesive.  Edges 
of  the  test  plate  were  smoothened  with  sandpaper  to  ensure  smooth  transition  of 
flow.  Silicon  sealant  was  applied  to  prevent  leakages.  This  window  section  would 
allow  flexing  of  test  plate  and  therefore  the  measurement  of  strain  of  the  plate 
during  the  experiment,  as  seen  in  Figure  35. 


Figure  35.  Test  Section 


Finally,  all  the  sections  were  assembled  using  clamps  and  sealant  to  form 
the  entire  duct,  3.05  m  in  length,  0.254  m  in  width  and  0.0318  m  in  height,  see 
Figure  36. 


Figure  36.  Entire  Duct 


2.  Flow  Circuit 

Due  to  the  requirement  to  create  a  laminar  flow  with  high  viscous  fluid,  a 
controllable  and  accurate  flow  circuit  was  constructed,  see  Figure  37. 
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Figure  37.  Flow  Circuit 


The  flow  circuit  was  powered  by  an  Baldor  VL3515  AC  motor  single  phase 
pump,  with  a  maximum  flow  rate  of  70  GPM.  Flow  diffusers  were  used  to  direct 
excessive  flow  away  from  the  duct,  see  Figure  38. 


Figure  38.  AC  Motor  Single  Phase  Pump 
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3. 


Instrumentation 


a.  Strain  Gages 


Vishay  Micro-Measurements  CEA-00-250UR-350  rosette  strain 
gages  were  used  to  measure  the  strain.  As  shown  in  Figure  39,  two  gages  were 
fitted  to  the  underside  of  the  test  plate. 


Figure  39.  Rosette  strain  gages 


b.  Data  Analyzer 

Measurements  from  the  rosette  strain  gages  were  acquired  by  a 
data  acquisition  system.  The  system  consisted  of  a  Pentium™  4,  2.4  GFIz,  512- 
MB  RAM  system  and  National  Instruments™  simultaneous  sampling 
multifunction  DAQ.  This  system  is  shown  in  Figure  40. 
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Figure  40.  Data  Acquisition  System 

The  system  had  a  16  bit  analog-to-digital  conversion  resolution  and 
was  capable  of  reading  a  total  of  16  channels  at  a  throughput  rate  of  up  to  250 
kS/s  per  channel,  which  was  appropriate  for  the  rate  of  testing  used  in  this  study. 
The  data-acquisition  process  was  controlled  using  the  NI-DAQmx  driver  software 
and  LabVIEW™  interactive  data-logging  software. 
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III.  FLUID  STRUCTURE  INTERACTION  MODELING 


A.  OVERVIEW 

In  order  to  understand  transient  response  of  the  test  plate  under  FSI 
conditions,  computational  studies  were  carried  out  in  ANSYS  13.0  Multi-field 
(MFX)  simulation.  Structural  analysis  of  the  test  plate  was  set  up  in  Transient 
Structural  analysis  system  while  the  flow  analysis  concerning  the  test  fluid  was 
set  up  in  Fluid  Flow  (CFX)  analysis  system.  Coupling  between  the  two  analysis 
systems  was  required  throughout  the  solution  to  model  the  FSI  between  test 
plate  and  fluid  as  time  progresses.  The  geometry  used  in  the  simulation 
consisted  of  the  fluid  body  and  the  test  plate  body  that  followed  the  dimensions  of 
the  actual  experimental  setup  (see  Figure  41). 


Figure  41.  Geometry  Used  in  FSI  Modeling 

Different  sizes  of  meshes  were  used  for  the  fluid  body  and  test  plate  body 
since  the  fluid  body  was  much  larger  than  the  test  plate  body.  The  program  was 
able  to  handle  the  mismatch  of  meshes  through  an  interpolation  algorithm.  In  all 
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cases,  the  test  plate  body  had  about  1,250  tetrahedron  elements  while  the  fluid 
body  had  about  500,000  tetrahedron  elements.  Boundary  conditions  applied  to 
the  test  plate  body  include  fixed  support  at  the  edges  and  fluid-solid  interface 
which  defined  the  interface  between  the  fluid  body  and  the  test  plate  body.  For 
the  fluid  body,  inlet  velocity,  outlet  pressure  and  wall  boundary  conditions  at  the 
interface  were  applied.  The  boundary  condition  at  the  interface  allowed  nodal 
velocity  values  to  be  received  from  the  Transient  Structural  solver  and  the 
sending  of  force  values  to  the  Transient  Structural  solver.  The  simulations  were 
conducted  for  a  duration  of  0.6s  and  with  a  time  step  of  0.005s.  Parameters  such 
as  nodal  displacement,  fluid  pressure,  structural  acceleration  and  velocity  were 
gathered  for  each  simulation.  The  studies  included  comparisons  of  results  of 
different  velocities,  densities,  and  Young’s  Moduli.  Propylene  Glycol  was  used  as 
the  test  fluid  in  all  the  cases.  The  baseline  case  was  set  up  as  close  as  possible 
to  the  actual  experiment.  Thus,  Al  1100  properties  were  used  and  a  uniform 
velocity  of  1  m/s  was  prescribed  at  the  inlet. 

B.  VARIATIONS  IN  INPUT  VELOCITY 

In  the  first  example,  the  inlet  velocity  was  varied  for  comparison.  Inlet 
velocities  of  0.5  m/s  and  2  m/s  were  modeled  and  compared  to  the  baseline 
case.  Results  for  nodal  displacement,  fluid  pressure,  structural  acceleration  and 
velocity  were  illustrated  in  Figures  42  to  45.  Nodal  displacements  for  both  fluid 
and  test  plate  were  similar,  thus  only  one  of  them  is  shown  here.  Undefined  data 
was  present  in  the  plot  in  Figure  43,  hence  producing  a  curve  that  was  not 
sinusoidal  in  nature. 
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Figure  42.  Nodal  Displacement  Z  for  Fluid  Body  for 
Different  Inlet  Velocities 


Figure  43.  Acceleration  Z  for  Test  Plate  for  Different  Inlet  Velocities 
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Figure  44.  Pressure  for  Fluid  Body  for  Different  Inlet  Velocities 
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Figure  45.  Velocity  W  for  Fluid  Body  for  Different  Inlet  Velocities 
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As  shown  in  the  preceding  group  of  figures,  in  all  three  cases,  vibration  of 
the  test  plate  damped  out  as  time  progressed  because  of  the  FSI  effect. 
Damping  out  started  at  about  0.25s  for  inlet  velocity  of  0.5  m/s,  0.30s  for  1  m/s 
and  beyond  0.60s  for  2  m/s.  It  was  observed  that  the  magnitude  of  displacement, 
acceleration,  pressure  and  velocity  at  the  interface  increased  with  inlet  velocity. 

C.  VARIATIONS  IN  DENSITY 

In  this  example,  the  density  of  the  test  plate  was  varied  for  the  simulations, 
while  all  other  material  properties  and  experimental  conditions  remained. 
Densities  of  1335  kg/m3  (half  of  Al  1100’s  density)  and  5420  kg/m3  (twice  of  Al 
1100’s  density)  were  modeled.  Results  for  nodal  displacement,  fluid  pressure, 
structural  acceleration  and  velocity  were  illustrated  in  Figures  46  to  49.  Nodal 
displacements  for  both  fluid  and  test  plate  were  similar,  thus  only  one  of  them  is 
shown  here.  Undefined  data  was  present  in  the  plot  in  Figure  47,  hence 
producing  a  curve  that  was  not  sinusoidal  in  nature. 
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Figure  46.  Nodal  Displacement  Z  for  Fluid  Body  for 

Different  Densities 
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Figure  47.  Acceleration  Z  for  Test  Plate  for  Different  Densities 
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Figure  48.  Pressure  for  Fluid  Body  for  Different  Densities 
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Figure  49.  Velocity  W  for  Fluid  Body  for  Different  Densities 


It  could  be  seen  from  the  results  that  as  density  increased,  the  FSI  effect 
decreased.  However,  for  all  the  cases,  differences  in  magnitude  of  the 
parameters  measured  were  not  significant.  The  vibration  for  all  cases  damped 
out  at  about  0.25s. 

D.  VARIATIONS  IN  YOUNG’S  MODULUS 

Next,  the  Young’s  Modulus  of  the  test  plate  was  varied  to  investigate  its 
effect  on  FSI,  while  all  other  material  properties  and  experimental  conditions 
remained  constant.  Young’s  Moduli  of  34.5  GPa  (half  of  Al  1100’s  Young’s 
Modulus)  and  134.9  GPa  (twice  of  Al  1100’s  Young’s  Modulus)  were  modeled  in 
the  simulations.  Results  for  total  mesh  displacement,  fluid  pressure,  structural 
acceleration  and  velocity  were  illustrated  in  Figures  50  to  53.  Nodal 
displacements  for  both  fluid  and  test  plate  were  similar,  thus  only  one  of  them  is 
shown  here.  Undefined  data  was  present  in  the  plot  in  Figure  51,  hence 
producing  a  curve  that  was  not  sinusoidal  in  nature. 
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Figure  50.  Nodal  Displacement  Z  for  Fluid  Body  for 
Different  Young’s  Moduli 


Figure  51 .  Acceleration  Z  for  Test  Plate  for  Different 

Young’s  Moduli 
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Figure  52.  Pressure  for  Fluid  Body  for  Different  Young’s  Moduli 


Figure  53.  Velocity  W  for  Fluid  Body  for  Different  Young’s  Moduli 
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As  illustrated  in  the  preceding  group  of  figures,  under  the  same  input 
velocity  for  fluid  body,  the  FSI  effect  was  greater  for  lower  Young’s  Moduli.  This 
was  exemplified  by  the  results  for  Young’s  Modulus  of  34.5  GPa.  The  magnitude 
of  parameters  measured  was  much  larger  than  the  other  two  simulations.  The 
FSI  effect  damped  out  at  about  0.25s  for  both  Young’s  Modulus  of  68.9  GPa  and 
134.9  GPa.  On  the  other  hand,  for  the  lowest  Young  Modulus  of  34.5  GPa,  the 
FSI  effect  did  not  showed  any  signs  of  damping  out  for  the  duration  of  simulation. 
These  results  were  consistent  with  the  physical  meaning  of  Young’s  Modulus  as 
a  measure  of  stiffness. 
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IV.  EXPERIMENTAL  RESULTS 


A.  OVERVIEW 

The  experiment  was  performed  using  the  designed  and  constructed  duct 
and  flow  circuit.  Many  pre-tests  were  conducted  to  detect  and  prevent  leakages 
and  to  ensure  that  the  flow  circuit  was  in  appropriate  working  conditions  for  the 
experiment.  Measurements  of  strains  were  recorded  for  a  period  of  120  sec  with 
a  sampling  rate  of  100/s  and  resulting  in  12,000  data  per  strain  gage.  The  gages 
and  pump  were  initiated  together  to  synchronize  the  measurements.  The  flow 
was  able  to  fill  up  the  whole  duct  and  became  steady  after  about  25  sec.  The 
pump  was  switched  off  after  the  flow  had  reached  an  equilibrium  state.  Due  to 
the  lack  of  instrumentation,  pressure  and  velocity  of  flow  were  not  measured. 
Five  tests  were  conducted  to  confirm  repeatability  of  the  testing  data. 

B.  RESULTS 

The  strain  values  were  measured  from  the  gages  on  the  underside  of  the 
test  plate,  with  the  positions  and  naming  shown  in  Figure  54. 


Figure  54.  Positions  and  Naming  of  Strain  Values 
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1.  Experimental  Results 

Figures  55  to  61  illustrated  a  representative  sample  of  the  experimental 
results  from  the  strain  gages’  measurement  for  one  of  the  runs.  Results  were 
consistent  for  all  the  runs. 


Figure  55.  £A  Response 


Figure  56.  £B  Response 
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Figure  57.  £c  Response 
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Figure  58.  £D  Response 
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Figure  59.  £e  Response 
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Figure  60.  £f  Response 
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Figure  61 .  Combining  Results  for  Strains 


From  the  preceding  group  of  figures,  the  following  trend  could  be 
observed  generally.  At  about  15  sec.,  the  strain  values  £a,  and  Ed  started  to 
increase  sharply  while  strain  values  £c  and  £f  decreased  sharply  as  the  fluid 
started  to  fill  up  the  duct.  Between  20  sec  to  60  sec,  the  strain  values  were 
oscillating  about  some  equilibrium  point  as  the  flow  approaches  steady  state.  At 
about  60  sec,  the  pump  was  switched  off  intentionally  to  observe  changes  in  the 
results.  The  strain  values  decreased  sharply  before  damping  out  as  the  flow 
slowed  down  from  60  sec.  onwards.  The  decay  of  £□  and  £f  were  the  greatest. 

The  strain  values  (£a  and  £c)  from  the  gage  at  center  of  the  plate  were 
smaller  than  the  strain  values  (£□  and  £f)  from  the  gage  nearer  to  the  boundary  in 
terms  of  magnitude,  as  expected.  Comparing  the  strains  along  x-axis  and  y-axis, 
the  magnitudes  at  each  gage  were  different  in  signs.  This  was  due  to  the  uneven 
surface  of  the  test  plate  where  local  warping  and  kinking  existed. 
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Lastly,  results  across  the  five  runs  were  consistent,  illustrating  the 
repeatability  of  the  experiment. 


2.  Comparison  to  Computational  Results 

Strain  data  from  the  previous  FSI  computational  modeling  of  the  baseline 
case  are  shown  from  Figures  62  to  65  for  comparison. 
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Figure  62.  £A  Response  from  Computational  Modeling  of  Baseline  Case 
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Figure  63.  £c  Response  from  Computational  Modeling  of  Baseline  Case 


Baseline 


Figure  64.  £p  Response  from  Computational  Modeling  of  Baseline  Case 
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Figure  65.  £f  Response  from  Computational  Modeling  of  Baseline  Case 

From  the  comparison  between  the  experimental  strain  values  and  strain 
values  from  the  computational  modeling  of  the  baseline  case,  it  was  observed 
that  the  order  of  magnitude  for  the  strain  values  differed  by  about  IxlO-1. 
Because  pressures  and/or  inlet  and  outlet  velocities  could  not  be  measured, 
more  quantitative  comparison  could  not  be  made.  In  addition,  for  £c  and  £f,  the 
trend  of  response  from  the  experiment  was  reasonably  similar  as  compared  to 
the  computational  values. 

To  explore  the  possibility  of  creating  larger  deflections  on  the  test  plate, 
the  current  model  was  modified  and  simulated  in  ANSYS  for  analysis.  The  fluid 
body  was  lengthened  for  in  the  -Z  direction  above  the  test  plate  in  an  attempt  to 
create  a  swirling  motion  in  a  cavity,  as  shown  in  Figure  66.  Two  models  with 
extensions  of  0.102m  and  0.203m  were  simulated,  respectively.  A  sample 
velocity  vectors  plot  at  0.03  sec  was  illustrated  in  Figure  67. 
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Figure  66.  Modified  Model  for  Simulation 


Figure  67.  Velocity  Vector  Plot  of  Modified  Model 

As  shown  in  Figure  67,  swirling  motion  did  occur  in  the  fluid  body  inside 
the  cavity  box.  Flowever,  as  illustrated  in  Figures  68  to  71,  an  examination  of  the 
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strain  values  at  the  test  plate  did  not  show  significant  increases  in  magnitude. 
Further  investigations  would  be  necessary  to  verify  these  results. 
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Figure  68.  Comparing  £a  Responses 
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Figure  69.  Comparing  £c  Responses 
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Figure  70.  Comparing  £d  Responses 
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Figure  71 .  Comparing  £f  Responses 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  study,  a  new  experimental  facility  was  designed  and  developed  for 
the  investigation  of  FSI  effect  between  a  metal  plate  and  fully  developed  laminar 
flow  through  a  duct.  Appropriate  dimensions  of  the  duct  were  obtained  through  a 
series  of  modeling  and  simulations.  The  construction  of  the  facility  was  carried 
out  after  a  careful  selection  of  materials.  The  experiment  was  performed  and 
strain  data  were  collected.  As  a  parametric  study,  fully  coupled  computational 
simulations  were  conducted  to  understand  the  transient  response  of  the  metal 
plate  under  FSI  conditions.  From  the  study,  it  was  found  that  the  magnitude  of 
displacement,  acceleration,  pressure  and  velocity  at  the  interface  increased  with 
inlet  velocity.  FSI  effect  was  also  observed  to  decrease  as  density  and  Young’s 
Modulus  increased. 

Results  from  the  experiments  proved  that  the  facility  developed  was 
capable  of  maintaining  constant  flow  through  the  duct  when  operating  with  a  fluid 
of  high  viscosity  and  high  reproducibility  of  results. 

An  additional  computational  study  was  conducted  to  examine  the 
likelihood  of  generating  larger  deflections  by  including  a  cavity  box  above  the  test 
plate.  Results  of  the  study  showed  that  a  swirling  motion  occurred  in  the  fluid 
body  in  the  cavity  box  as  flow  past  through,  but  it  did  not  increase  deflection  of 
the  test  plate.  This  could  be  verified  in  further  investigations. 

Other  future  studies  could  include  more  measurement  techniques  to 
investigate  other  critical  FSI  parameters  such  as  displacement,  velocity,  fluid 
pressure  and  frequency.  For  example,  using  non-contact  inductive  distance 
sensors  or  digital  correlation  technique  for  displacement  measurement  and 
differential  pressure  sensors  for  the  assessment  of  fluid  pressure  and  velocity. 
Investigations  with  other  plate  materials,  different  plate,  other  duct  geometries 
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and  fluids  with  different  viscosities  could  be  conducted  to  gain  more  insight  into 
the  FSI  effect,  hence  building  a  database  of  results  that  could  be  easily  validated 
and  verified. 
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